Separation Of Variables Research Articles

Signature of f(R) gravity via Lemaître–Tolman–Bondi inhomogeneous perturbations

We analyze inhomogeneous cosmological models in the local Universe, described by the Lemaître–Tolman–Bondi (LTB) metric and developed using linear perturbation theory on a homogeneous and isotropic Universe background. Focusing on the different evolution of spherical symmetric inhomogeneities, we compare the Λ\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\Lambda $$\\end{document}LTB model, in which the cosmological constant Λ\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\Lambda $$\\end{document} is included in the LTB formalism, with inhomogeneous cosmological models based on fR\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$f\\left( R\\right) $$\\end{document} modified gravity theories viewed in the Jordan frame. We solve the system of field equations for both inhomogeneous cosmological models adopting the method of separation of variables: we integrate analytically the radial profiles of local perturbations, while their time evolution requires a numerical approach. The main result of the analysis concerns the different radial profiles of local inhomogeneities due to the presence of a non-minimally coupled scalar field in the Jordan frame of fR\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$f\\left( R\\right) $$\\end{document} gravity. While radial perturbations follow a power-law in the Λ\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\Lambda $$\\end{document}LTB model, Yukawa-like contributions appear in the fR\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$f\\left( R\\right) $$\\end{document} theory. Interestingly, this latter peculiar behavior of radial profile is not affected by the choice of the fR\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$f\\left( R\\right) $$\\end{document} functional form. The numerical solution of time-dependent perturbations exhibits a non-diverging profile. This work suggests that investigations about local inhomogeneities in the late Universe may allow us to discriminate if the present cosmic acceleration is caused by a cosmological constant term or a modified gravity effect.

Integrability of Dirac equations in static spherical space-times

We consider the Dirac equations in static spherically-symmetric space-times, and we present a type of spinor field whose structure allows the separation of elevation angle and radial coordinate in very general situations. We demonstrate that after such a separation of variables the Dirac equations reduce to two equations that can always be integrated, at least in principle. To prove that ours is a fully-working method, we find an explicit exact solution in the special case of the de Sitter universe.

A novel approach to model nance drying considering its shrinkage

Air drying is a viable preservation method to boost the accessibility of nance, a seasonal drupe with very short shelf life; however, information about nance drying is scarce, especially from a mass transfer (MT) standpoint. Besides, the lack of simple models considering the stone fraction (radius of inner stone/radius of nance) and the mesocarp shrinkage in drupes has caused most investigations rely on sphere and cube geometries, contributing to the unreliable estimation of MT properties. This study developed a discrete state-space (SS) formulation with analytical solution to describe nance drying which, unlike a variable separation (VS) solution, can consider a variable stone fraction, caused by outer boundary moving while flesh-pit interface remains still. Nance fruits were air dried (60, 70 and 80 °C) to obtain their moisture kinetics and shrinking behavior. Drying data were used to estimate water diffusivity (D) under rigid (VS and SS solutions) and shrinkable (SS and numerical solution) mesocarp assumptions. Nances lost about 60% of their initial volume at the end of drying, while stone fraction increased from 0.44 to 0.62. D values ranged from 5.3 to 6.7 (×10–10) m2/s with SS approach and variable stone fraction. The major source of error in D was not the formulation (SS or VS) when a rigid mesocarp was considered, but the values assigned to the nance radius (length for diffusion) and stone fraction, an initial or average value along process. SS diffusivities did not differ to those estimated with the full numerical solution under variable stone fraction (p>0.05), while VS parameters were overestimated by 44–46 % as a minimum. SS formulation can be applied to hollow cylinders and other geometries with 2D or 3D MT, while able to include non-isotropic shrinkage and evolving external conditions, a problem normally reserved to pure numerical implementations, offering substantial advantages over VS equations.

A Study of Singular Similarity Solutions to Laplace’s Equation with Dirichlet Boundary Conditions

A Study of Singular Similarity Solutions to Laplace’s Equation with Dirichlet Boundary Conditions

An analytical fast field solver for surface permanent magnet co‐axial magnetic gears

AbstractA mathematical model of co‐axial magnetic gears is given for the prediction of magnetic field distribution and the magnetic forces. The modelling approach is based on the separation of variables method (SVM), Sturm–Liouville problem, and spectral analysis. This method is computationally more efficient and has the ability to account for finite permeability of magnetic saliency. An finite element method (FEM) validation is performed for a case study, and the results are observed to be almost matched. The computational time of the analytical model has been almost in the same order of that of the FEM for a specified rotor position (both in a matter of a second or less); however, when 100 number of rotor positions were accounted in the simulation, the computational time of FEM was increased to 82 s, while the computational complexity of SVM was almost unaffected.

Relativistic Massive and Massless Scalar Fields Bound to The Bumblebee Gravity's Black Hole

In this letter, we examine massive and massless scalar quasibound states bound in the static spherically symmetric black hole solution of the Einstein-Bumblebee space-time. We start with constructing the governing relativistic fields' dynamical equation, i.e. the Klein-Gordon equation, component by component. With the help of the ansatz of separation of variables, we find the exact solution of the angular part in terms of Spherical Harmonics while the radial exact solution is discovered in terms of the Confluent Heun function. The quantization of the quasibound state is done by applying the polynomial condition of the Confluent Heun function that reveals a complex-valued energy levels expression for the case of massive scalar, where the real part is the scalar particle's energy and the imaginary part represents the quasibound state's decay. And for a massless scalar, a pure imaginary energy levels is obtained. The quasibound states, thus, describe decaying relativistic states bound in the black hole's gravitational potential well. At the end, we investigate the Hawking radiation of the static Bumblebee black hole's horizon by deriving and investigating the radiation distribution function.

Analysis of Hydrodynamic Loads by a Vertical Hollow Cylindrical Structure in A Two-Layer Fluid of Finite Depth

This paper aims to study the effect of surge radiation by a vertically hollow cylinder floating in a two-layer fluid of finite depth. Since most of the important characteristics of the oscillating water column are preserved by the hollow cylindrical structure, so the proposed device can be considered as an oscillating water column (OWC) which is one of the important wave energy device. In this two-layer fluid system, the upper layer is bounded above by a free surface and the lower layer is bounded below by a solid impermeable bottom, and the submerged hollow cylindrical structure allows it to oscillate in a surge mode of motion in the upper layer. On the assumption that wave amplitudes are small in comparison to wavelengths, hence the linear water wave theory is applied to formulate boundary value problems by dividing whole domain into 2 sub-domains. The formulated boundary value problems for surge radiated potentials in each sub-domain; we execute to obtain solutions in each sub-domain by using variable separation and matched eigenfunction expansion methods. With the help of obtained surge radiated potentials, we derive the radiation forces by the device in terms of added mass and damping coefficients. A set of influences of different parameters of the device on the added mass and damping coefficients are demonstrated in both surface and internal wave modes of motion. Hydrodynamic coefficients oscillate highly near a particular frequency, which leads to the resonance phenomena. In addition to that, the 3D plots of the free surface elevation in surface and internal wave modes are presented. A density ratio γ = 0.97 has a higher oscillation than a density ratio γ = 0.93, 0.95. High oscillations occur around a particular frequency ω = 3.83 due to the resonance phenomenon in which the waves incoming match with the natural motion of the object. HIGHLIGHTS A floating hollow vertical cylindrical structure is considered in a two-layer fluid system having finite depth The separation of variables and matched eigenfunction expansion methods are utilized to obtain the solutions to the boundary-value problems The diffracted velocity potentials under the surge mode of motion are evaluated in the identified regions The added mass, damping coefficients, free surface and interface elevations in surface and internal wave modes are evaluated with different set of parameters of the device A number of results are presented regarding the effect of parameters of the device on the added mass, damping coefficients, free surface and interface elevations It is observed that the density ratio of the two-layer fluid has significant effect on the added mass, damping coefficient and the elevations GRAPHICAL ABSTRACT

New Multilinear Variable Separation Solutions of the (3 + 1)-Dimensional Burgers Hierarchy

The multilinear variable separation (MLVS) approach has been proven to be very useful in solving (2 + 1)-dimensional integrable systems. Taking the (3 + 1)-dimensional Burgers hierarchy as an example, we extend the MLVS approach to a whole family of (3 + 1)-dimensional Burgers hierarchy. New exact solutions and universal formulas are obtained, which lead to abundant (3 + 1)-dimensional coherent structures. In particular, two ring-type soliton molecules and their interactions are shown in detail. We also generalize the MLVS results of the (3 + 1)-dimensional Jimbo–Miwa (JM) equation and modified JM equation.

Stationary coupled KdV systems and their Stäckel representations

AbstractIn this article, we investigate stationary coupled Korteweg–de Vries (cKdV) systems and prove that every ‐field stationary cKdV system can be written, after a careful reparameterization of jet variables, as a classical separable Stäckel system in different ways. For each of these parameterizations, we present an explicit map between the jet variables and the separation variables of the system. Finally, we show that each pair of Stäckel representations of the same stationary cKdV system, when considered in the phase space extended by Casimir variables, is connected by an appropriate finite‐dimensional Miura map, which leads to an ‐Hamiltonian formulation for the stationary cKdV system.

3D electromagnetic analytical model for radial-flux eddy-current couplers

Purpose This paper aims to develop a new 3D analytical model in cylindrical coordinates to study radial flux eddy current couplers (RFECC) while considering the magnetic edge and 3D curvature effects, and the field reaction due to the induced currents. Design/methodology/approach The analytical model is developed by combining two formulations. A magnetic scalar potential formulation in the air and the magnets regions and a current density formulation in the conductive region. The magnetic field and eddy currents expressions are obtained by solving the 3D Maxwell equations in 3D cylindrical coordinates with the variable separation method. The torque expression is derived from the field solution using the Maxwell stress tensor. In addition to 3D magnetic edge effects, the proposed model takes into account the reaction field effect due to the induced currents in the conducting part. To show the accuracy of the developed 3D analytical model, its results are compared to those from the 3D finite element simulation. Findings The obtained results prove the accuracy of the new developed 3D analytical model. The comparison of the 3D analytical model with the 2D simulation proves the strong magnetic edge effects impact (in the axial direction) in these devices which must be considered in the modelling. The new analytical model allows the magnetic edge effects consideration without any correction factor and also presents a good compromise between precision and computation time. Practical implications The proposed 3D analytical model presents a considerably reduced computation time compared to 3D finite element simulation which makes it efficient as an accurate design and optimization tool for radial flux eddy current devices. Originality/value A new analytical model in 3D cylindrical coordinates has been developed to find the electromagnetic torque in radial flux eddy current couplers. This model considers the magnetic edge effects, the 3D curvature effects and the field reaction (without correction factors) while improving the computation time.

Numerical modelling of thermal jet assisted rock cutting with double PDC cutters

AbstractPreconditioning of rock for drilling operations is a potential method to facilitate the mechanical breakage and mitigate the tool wear. This paper numerically investigates one such preconditioning technique, namely, the thermal jet assisted rock cutting. For this end, a numerical method for solving the governing thermo‐mechanical problem is developed and validated. The continuum approach is chosen to describe the rock failure, being based on a damage‐viscoplasticity model with the Drucker–Prager yield surface and the modified Rankine surface as the tensile cut‐off. In the damage part, separate scalar damage variables are employed for tension and compression. The cutters are modelled as rigid bodies and their interaction with the rock is modelled by imposing contact constraints by the forward increment Lagrange multiplier method, which is compatible with the chosen explicit time integration scheme. A damage‐based erosion criterion is applied to remove the contact nodes surrounded by heavily damaged elements. The thermal jets are modelled with a moving external heat flux boundary condition. The global thermo‐mechanical problem is solved with a staggered approach explicitly in time while applying mass scaling to increase the critical time step. The novel 3D numerical simulations involving two cutters demonstrate the capabilities of the method with a special emphasis on cutter‐thermal jet configurations. Therefore, the present method provides a potential tool for the bit design in thermal jet assisted drilling.

Separation of variables for scalar-valued polynomials in the non-stable range

Any complex-valued polynomial on (Rn)k decomposes into an algebraic combination of O(n)-invariant polynomials and harmonic polynomials. This decomposition, separation of variables, is granted to be unique if n≥2k−1. We prove that the condition n≥2k−1 is not only sufficient, but also necessary for uniqueness of the separation. Moreover, we describe the structure of non-uniqueness of the separation in the boundary cases when n=2k−2 and n=2k−3.Formally, we study the kernel of a multiplication map ϕ carrying out separation of variables. We devise a general algorithmic procedure for describing Ker ϕ in the restricted non-stable range k≤n<2k−1. In the full non-stable range n<2k−1, we give formulas for highest weights of generators of the kernel as well as formulas for its Hilbert series. Using the developed methods, we obtain a list of highest weight vectors generating Ker ϕ.

The Kerr–Bumblebee exact massive and massless scalar quasibound states and Hawking radiation

In this letter, we will focus on the Klein–Gordon equation with rotating axially symmetric black hole solution of the Einstein–Bumblebee theory, so called the Kerr–Bumblebee black hole, as its 3 + 1 background space-time. We start with constructing the covariant Klein–Gordon equation component by component and with the help of the ansatz of separation of variables, we successfully separate the polar part and found the exact solution in terms of Spheroidal Harmonics while the radial exact solution is discovered in terms of the Confluent Heun function. The quantization of the quasibound state is done by applying the polynomial condition of the Confluent Heun function that is resulted in a complex-valued energy levels expression for a massive scalar field, where the real part is the scalar particle’s energy while the imaginary part represents the quasibound stats’s decay. And for a massless scalar, a pure imaginary energy levels is obtained. The quasibound states, thus, describe the decaying nature of the relativistic scalar field bound in the curved Kerr–Bumblebee space-time. We also investigate the Hawking radiation of the Kerr–Bumblebee black hole’s apparent horizon via the Damour–Ruffini method by making use the obtained exact scalar’s wave functions. The radiation distribution function and the Hawking temperature are successfully obtained.

Control co-design optimization of nonlinear wave energy converters

This paper presents a study in which both the control and the shape of a Wave Energy Converter (WEC) are optimized simultaneously. A heaving point absorber WEC is assumed. To optimize the shape of the WEC’s buoy, nonlinear hydrodynamics need to be evaluated. One main contribution of this paper is the integration of nonlinear hydrodynamics and nonlinear control during the optimization of the WEC’s buoy shape. This approach is referred to as Control Co-Design (CCD). In this work, we present a control co-designed nonlinear heaving point absorber WEC that leverages the nonlinear dynamic, static Froude–Krylov (FK) forces to maximize power extraction. The nonlinear FK forces are approximated using a variation of the algebraic solution; the hydrodynamic forces of the body are computed using an analytic formulation leveraging the methods of eigenfunction expansion and separation of variables. The nonlinear geometry of the buoy is modeled as a series of inclined panels; the inclination angles are optimized to arrive at the optimal shape. The performance of the optimized shape is compared to that of a nonlinear spherical WEC. It is found that an average of 20% improvement is achieved by the optimized geometry over the spherical device.

On correlation functions for the open XXZ chain with non-longitudinal boundary fields: The case with a constraint

This paper is a continuation of [J. Phys. A: Math. Theor. 55, 405203 (2022)], in which a set of matrix elements of local operators was computed for the XXZ spin-1/2 open chain with a particular case of unparallel boundary fields. Here, we extend these results to the more general case in which both fields are non-longitudinal and related by one constraint, allowing for a partial description of the spectrum by usual Bethe equations. More precisely, the complete spectrum and eigenstates can be characterized within the Separation of Variables (SoV) framework. One uses here the fact that, under the constraint, a part of this SoV spectrum can be described via solutions of a usual, homogeneous, TQ-equation, with corresponding transfer matrix eigenstates coinciding with generalized Bethe states. We explain how to generically compute the action of a basis of local operators on such kind of states, and this under the most general boundary condition on the last site of the chain. As a result, we can compute the matrix elements of some of these basis elements in any eigenstate described by the homogenous TQ-equation. Assuming, following a conjecture of Nepomechie and Ravanini, that the ground state itself can be described in this framework, we obtain multiple integral representations for these matrix elements in the half-infinite chain limit, generalizing those previously obtained in the case of longitudinal boundary fields and in the case of the special boundary conditions considered in [J. Phys. A: Math. Theor. 55, 405203 (2022)].

A novel (2+1)-dimensional Sawada-Kotera type system: multisoliton solution and variable separation solution

A novel (2+1)-dimensional Sawada-Kotera type system: multisoliton solution and variable separation solution

A new analytical method for transient heat conduction in composite disks applied to thermal plug and abandonment of oil wells

A new hybrid method for transient heat conduction problems is developed and applied to simulate a Plug and Abandonment (P&A) operation of oil wells. An oil well is approximated by concentric disks in a one-dimensional configuration, which allows for use of polar coordinates. The application of the Separation of Variables Method (SVM) is used as the analytical framework for the solution of the conductive heat transfer arising from a volumetric heat source located in the center of the disks. The SVM is able to solve only time-independent boundary conditions. However, using the Duhamel's theorem, the solution determined with the SVM can be used to achieve the solution when both time-dependent boundary conditions and internal heat generation are prescribed. Lastly, for verification purposes, a commercial software that solves the transient temperature field by means of numerical procedures, providing reliability to the analytical method proposed in this work.

Soliton solutions for NLS model via the separation of variables method

Soliton solutions for NLS model via the separation of variables method

Hydrodynamic performance of vertical cylindrical wave energy absorbers in front of a vertical wall

This paper studies the hydrodynamic performance of vertical cylindrical absorbers in front of a vertical wall. All the absorbers are independent of each other and restricted to only the heave motion. Based on a linear potential flow theory, an analytical solution is developed for the problems of wave diffraction and radiation by absorbers. In the solving procedure, the hydrodynamic problem is first transformed into an equivalent problem in an open water domain using the image principle. The number of absorbers in the equivalent problem is twice that in the real problem, and the plane layout is symmetric about the original vertical wall. The velocity potential of the fluid domain is obtained using the method of variable separation, and the unknown expansion coefficients in the velocity potential are determined by the matched boundary conditions. The heave excitation force, added mass, radiation damping, motion response, and energy capture width of the absorbers are calculated. Case studies are presented to show the effects of the wall reflection and hydrodynamic interaction on the energy extraction performance of the wave energy converter (WEC) system. Subsequently, the WEC performance under the action of irregular waves is analyzed by considering an incident wave spectrum, and the mean annual absorbed power of the device is estimated by considering the wave data statistics at the actual sites. The results indicate that when the wave motion resonates with the absorber motion, the energy extraction performance of the absorbers is significantly improved. The performance of the absorbers can be effectively improved when the structures are close to the antinodes of a standing wave field. By designing a reasonable plane layout, the hydrodynamic interaction can play a constructive role in the performance of the WEC system.

A mathematical model of tissue axial and radial diffusion in the microvasculature for intravascular microscopy and phosphorescence quenching data

This study aims to extend earlier Krogh Cylinder Models of an oxygen profile by considering axial diffusion and analytically solving Fick's Law Partial Differential Equation with novel boundary conditions via the separation of variables. We next prospectively collected a total of 20 animals, which were randomly assigned to receive either fresh or two-week-old stored red blood cell (RBC) transfusions and PQM oxygen data were measured acutely (90 min) or chronically (24 hours). Transfusion effects were evaluated in vivo using intravital microscopy of the dorsal skinfold window chamber in Golden Syrian Hamsters. Hamsters were initially hemorrhaged by 50% of total blood volume and resuscitated one-hour post hemorrhage. PQM data were subsequently collected and fit the derived 2D Krogh cylinder model. Systemic hemodynamics (mean arterial pressure, heart rate) were similar in both pre and post-transfusion with either stored or fresh cells. Transfusion with stored cells was found to impair axial and radial oxygen gradients as quantified by our model and consistent with previous studies. Specifically, we observed a statistically significant decrease in the arteriolar tissue radial oxygen gradient after transfusion with stored RBCs at 24 hours compared with fresh RBCs (0.33 ± 0.17 mmHg μ m-1 vs, 0.14 ± 0.12 mmHg μ m-1; p = 0.0280). We also observed a deficit in the arteriolar tissue oxygen gradient (0.03 ± 0.01 mmHg μ m-1 fresh vs. 0.018 ± 0.007 mmHg μ m-1 stored; p = 0.0185). We successfully derived and validated an analytical 2D Krogh cylinder model in an animal model of microhemodynamic oxygen diffusion aberration secondary to storage lesions.